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^ ■ ABSTRACT 
O 

' We present a particle physics quintessence model which agrees well with existing cosmological 

data, including the position of the acoustic peaks. This model has an inverse power law potential 
2 ■ (IPL) with n = 18/7 ~ 2.57 and it gives tfe// = —0.75, an acoustic scale Ia = 307 and a density 

c/3 . contrast erg = 0.95. 

Models with n > 1 have been said to be disfavored by the analysis of the acoustic peaks. 
However, the results are not correct. The main reason is that the tracker approximation has 
j_j ■ been used in deriving the IPL constrains and for n < 5 the scalar field has not reached its tracker 

^ value by present day. 

The model can be derived from particle physics, using Affleck-Dine-Seiberg "ADS" superpoten- 
tial, for a non-abelian gauge group with Nc = 8,Nf = 1. The advantage of having Nf = 1 is 
that there is only one degree of freedom below the condensation scale given by the condensate 
(quintessence) cfP' =< QQ > field. The condensation scale is at IGeV a very interesting scale 
since it connects the quintessence "Q" with the standard model "SM" scale. The similarity in 
energy scales between Q and SM scale gives an "explanation" to the coincidence problem. The 
fact that only recently the universe is accelerating is a natural consequence of the Q scale and 
the evolution of (h. 



^e-mail: macorra@fisica.unani.nix 
^e-mail: gabriel@fis.unani.mx 



The Maxima and Boomerang observations on the cosmic radiation microwave background 
("CMBR") and the superonovae project SNla 0] have lead to conclude that the universe is flat 
and it is expanding with an accelerating velocity. These conclusions show that the universe is 
now dominated by an energy density r2<^o = 0.7 ± 0.1 (the subscript "o" refers to present day 
quantities) with negative pressure. The SNla data requires an equation of state w^o < —2/3 Q 
while recent analysis on the CMBR peaks constrains the models to have Weff = —0.82't'\\ ||5|, ^, 
where We// is an average equation of state. This energy is generically called the cosmological 
constant. Structure formation also favors a non- vanishing cosmological constant consistent with 
SNla and CMBR observations An interesting parameterization of this energy density is 
in terms of a scalar field with gravitationally interaction only called quintessence [10|. The 
evolution of scalar field has been widely studied and some general approach con be found in 
|l^ , |l5| , The evolution of the scalar field cf) depends on the functional form of its potential 
V{(j)) and a late time accelerating universe constrains the form of the potential [15|. 

One of the simplest and most interesting quintessence potentials are the inverse power law (IPL) 
|11|. In some special cases they can be derived from non-abelian gauge theories [12, 13, ^ and 
we can also have consistent models with a gauge coupling unified with the standard model (SM) 
couplings [13|. 

From the CMBR analysis it has been inferred that IPL with n < 1 are disfavored [^,0. However, 
in most cases one assumes a constant w given by the tracker value wtr = — 2/(2 + n) and for IPL 
models with n < 5 the tracker solution |1C] is not a good approximation to the numerical (exact) 
solution since the field has not reached its tracker value by present time. This fact implies that 
we cannot use the tracker solution (i.e. the constant wtr) for the evolution of (j) in determining 
the acoustic peaks. It is no surprise that the values of the acoustic peaks differ greatly if we use 
the tracker solution approximation or we evolve the quintessence field from its initial conditions. 
Models with 1 < n < 2.5 (for > 0.25) are therefore still phenomenologically viable. 

Tracker fields have the advantage that the value of w^o does not depend on the initial condition. 
In fact one can have more then 100 orders of magnitude on the initial conditions il^j ^3| and the 
value of Wtf)o will not change. In our case this is no longer so since the scalar field has not reached 
its tracker value and depends slightly on Q^^i, so there is a dependence on but there is 
no fine tuning required on the initial conditions (we do not need to adjust the initial condition 
more than one significant figure). Furthermore, in our model, derived from non-abelian gauge 
theories, we can determine the initial conditions in terms of the number of degrees of freedom of 
the system. So we do not need 100 orders of magnitude independence of initial conditions since 
they are well motivated and given within one order of magnitude at most. Of course, any model 
which requires a fine tuning on the initial conditions would not be theoretically acceptable but 
this is not our case. In both cases, tracker and our model, one still has the coincidence problem 
since the scale Ac has to be tuned so that ~ 0.7 with ho 0.7 today. 

There are two constrains on the equation of state parameter w^i,, one coming from direct obser- 
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vations SNla which sets un upper hmit, w^j^o < —2/3 Q and the other is indirect and comes from 
numerical analysis of the CMBR data and gives a smaller value Weff = — [^]. Notice 
that the CMBR data gives a more negative w^j, than the SNla one but it is an average equation 
of state (from last scattering to present day) while the SNla result gives an in recent times. 
In IPL with n < 5, where the quintessence field has not reached its tracker value yet, one has 
always W(f,o larger than Weff in good agreement with SNla and CMBR data. For IPL models 
it was shown |18] that depends on n and the initial condition il^j. If we want W(j,o < — 2/3 



assuming an > 0.25 IPL models require an n to be less 2.74 |18] assuming no contribution 
from radiation at present time. If we include radiation with Clro = 4.17 x 10~^hg then the value 
of n will decrease slightly, e.g for Q^pi = 0.25 we have Wt^o < —2/3 for n < 2.5. Larger values 
of flffii allow larger values of n, however we would not expect to have much larger since a 
"reasonable" amount of energy must go into the standard model of elementary particles "SM". 
These results set an upper value of n but there is still room for models with 1 < n < 2.5 and if 
we take Q^j, > 0.3 then the value of n can be as large as n < 2.66. 

The CMBR constrain on can be studied from the position of the third acoustic peak. The 
position of the third CMBR peak has been found to be not very sensitive to the different 
cosmological parameters and it is a good quantity to obtain the acoustic scale I a [§]• The 
acoustic scale Ia, which sets the scale of the peaks, derived from the third acoustic peak is 

/a = 316 ±8 (1) 

were we have taken ^3 = 845^25 Q (see below for the definition of Ia)- 

Here we will present a model with the largest value of n that is still in agreement with the 
observational cosmological data Q and that can be nicely derived from particle physics. Since, 
for a non-abelian gauge group (see below) we have n = 2 + 4Nf/(Nc — Nf) = 2 + 4/{z — 1), 
where Nf is the number of fiavors of the gauge group SU{Nc) and z = Nc/Nf, we can see 
that n decreases with increasing z. The requirement on n to be smaller than 2.66 (2.5) for 
^if>i > 0.3 (0.25), in order to give the correct phenomenology and the observed values of the 
acoustic scale and present day tw^o, implies that 2; > 9 (7) or z Nj < Nc- As we will see later 
the number of condensates of the gauge group SU{Nc) is given in terms of Nf. Therefore, the 
model with the least number of condensates has Nf = 1 and the smallest gauge group would be 
Nc = 8. Furthermore, in string compactification for the heterotic string, the gauge group has at 
most rank 8 (as SU{8)). The value of n for SU{8) with AT^ = 1 is n = 18/7 = 2.57. This model 
represents the limiting acceptable model. 

The acoustic scale gives Ia = 307 and Wcpo = —0.68, which gives a good prediction of the 
acoustic peaks within observational limits. This acoustic scale should be compared with the 
tracker solution lAtr = 281 and wtr = —2/(2 + n) = —0.44 and the cosmological constant 
lACte = 315 (i.e. wcte = — 1) for the same initial conditions. We see that the tracking solution 
is not a good approximation since w^po differs by more than 38% and the acoustic scale Ia by 9% 
from the numerical solution of the scalar field, discrepancy large enough to rule out the model. 
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Since our model has n < 5 the quintessence field has not reached its tracker value yet and the 
coincidence problem is not solved. However, there is a clear connection between the model con- 
densation scale Ac = IGeV and the standard model scale. So, we could think of the "solution" 
to the coincidence problem as the following: The scale of quintessence should not be given by 
today's energy density but by the condensation scale Ac. The natural value of this scale is that 
of the standard model. The subsequent evolution of the quintessence field is determined by 
the solution of the Friedmann-Robertson- Walker equations and the fact that only recently the 
universe is accelerating is a natural consequence of the quintessence dynamics starting at Ac. 

The model with n = 18/7 can be easily obtained from a SU{8) non-abelian gauge group with 
Nf = 1 number of (chiral + antichiral) fields in the fundamental representation and with a 
condensation scale Ac = IGeV, quite an interesting scale. Above the condensation scale the 
gauge coupling constant is small and the elementary fields are massless. At the condensation 
scale there is a phase transition, the gauge coupling constant becomes strong, and binds the 
elementary fields together forming meson fields. In this model there is only one degree of freedom, 
(j), in the confined phase and the Affeck-Dine-Seiberg "ADS" superpotential ||l^ obtained is 
therefore exact. 

At the beginning we have particles of the standard model (SM) and the quintessence model 
(Q). All fields, SM and Q model, are massless and redshift as radiation until we reach the 
condensation scale Ac of Q group. Below this scale the fields of the quintessence gauge group 
will dynamically condense and we use ADS potential to study its cosmological evolution. The 
ADS potential is non-perturbative and exact (it receives no quantum corrections) and it is 
given for a non-abelian SU{Nc) gauge group with Nf (chiral + antichiral) massless matter fields 
by 

W = {N,- Nf){ )yiN.-Nr) (2) 

det < QQ > 

where bo = 3Nc — Nf is the one-loop beta function coefficient. The scalar potential in global 
supersymmetry is V = \ W^\'^, with = dW/d(j), giving |12] 



V = c2A^+"0-'^ (3) 

where we have taken det < QQ >= li-'^Ii'p'j, c = 2Nf, n = 2 + 4jy^yj and Ac is the 
condensation scale of the gauge group SU{Nc). Our model has Nf = 1, Nc = 8 and n = 18/7. 
There are no baryons since Nf < Nc and there is only one degree of freedom below Ac which is 
the condensate cj) =< QQ >. 

The cosmological evolution of (j) with an arbitrary potential V{(j)) can be determined from a 
system of differential equations describing a spatially flat Friedmann-Robertson- Walker universe 
in the presence of a barotropic fluid energy density p-y that can be either radiation or matter, 
are 

H = -^{pj+Py + 'j)^), 
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p = -3H{p + p), (4) 

where H is the Hubble parameter {H = WOh km Mpc~^s~^), f = df/dt, p (p) is the total 
energy density (pressure) and we are setting the reduced Planck mass = l/SvrG = 1. 

Solving eqs.0 we have that the energy density of the Q group fi,^ drops quickly, independently 
of its initial conditions, and it is close to zero for a long period of time, which includes nucle- 
osynthesis (NS) if Ac is larger than the NS energy A^s (or temperature Tjys = 0.1 — lOMeV), 
and becomes relevant only until very recently On the other hand, if Ac < A^vs then the NS 
bounds on relativistic degrees of freedom must be imposed on the models. Finally, the energy 
density of Q grows and it dominates at present time the total energy density with the il^o ~ 0.7 
and a negative pressure w^o < — 2/3 leading to an accelerating universe 

The value of the condensation scale in terms of Hq is ||l^, |l8[ 

''=[-^) 

The approximated value can be obtained since one expects, in general, to have y^^o ^ 1 for 
a model with $7,^0 = 0.7 and w^po < —2/3. The magnitude order of the condensation scale is 
therefore A, = H^o^ '•^+''1 

In our model, the cosmological evolution requires a condensation scale Ac = IGeV in order to 
give = 0.7 ± 0.1 with a Hubble parameter ho = .65 ± .7 at present time. Since Ac » A^vs 
the energy density at NS is 0,(j){NS) <C 1 and there is no constrain from nucleosynthesis on the 
model. 

In order to set the initial conditions for (j) we will assume that all relativistic degree of freedom 
(MSSM and Q) had the same fraction of energy density at high energies, when all fields were 
massless. The initial conditions could be set at the unification scale or at the reheating tem- 
perature Trh. If Trh is larger then the supersymmetric masses (i.e. Trh > W^GeV), which is a 
natural assumption, then all degrees of freedom (MSSM and Q) would be relativistic at T^h and 
each degree of freedom would have the same energy density (assuming the standard reheating 
process which is gauge blind). Therefore, the initial energy density conditions would be exactly 
the same at the unification scale or reheating temperature. 

The MSSM has gsmi = 228.75 while the Q group has gqi = (l+7/8)(2(iV2 -l)+2iVjiVc) = 266.25 
de grees of freedom, where ga — YjuBosons + 7 /SYj^FeTniions. Taking into account that some 
fields become massive at lower energies, we can determine the energy density at an arbitrary 



energy scale A and it is given by [|13| 



9smf + 9Qf{9smf9Qi/9smi9QfY^^ 
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where (7smji 9smf, 9Qi, 5q/ are the initial (i.e. at high energy scale) and final (i.e. at A) standard 
model and Q model relativistic degrees of freedom, respectively. Taking A = Ac = IGeV the 
MSSM has gsmf = 10.75 and if the Q group is still super symmetric at Ac it has gqj = gqi 
and O0j(Ac) = 0.29. If Q is no longer supersymmetric gqf = gqij'^ and r2(^j(Ac) = 0.34. So a 
reasonable choice for the initial conditions is r20j(Ac) = 0.3. 

We show in fig.(||) the evolution of and as a function of = Log{a), with a the scale factor, 
for initial conditions Jl^, = 0.5, 0.3, 0.2 short-dashed, long-dashed and solid lines, respectively. 
We see that if^o decreases for larger initial condition Q^f^i. For il^j = 0.9, 0.5, 0.3, .2 one finds 
w^o = —0.95, —0.8, —0.68, —0.61 and an acoustic scale Ia = 314, 313, 307, 303, respectively. It 
is no surprise that for = 0.2 the value of = —0.61 lies outside the observed range 
W(f,o < —2/3. This is because the model we are working with (i.e. n=18/7=2.57) gives almost 
the limiting value of W(i,o = —2/3 for Jl^j = 0.3 ( W(i,o increases for smaller ^l^jji) and that was 
the reason for using this model. However, as long as we take > 0.27 the model satisfies the 
cosmological constrains. For scalar fields that have reached its tracker value by present day (i.e. 



71 > 5) the initial value of fi^j is not constrained |1C] since it can vary for more than 100 orders 
of magnitude and the value of w^o will be the same 'u^f^o — — ^/(^ ~^ ^) [|10|1 , however for ti > 5 
one has > —0.28 which is too large. 

In previous works |1^, 18 1 we have studied quintessence models that are unified with the standard 



model gauge groups, i.e. the gauge coupling constant of all gauge groups is the same at the 
unification scale. In the model we are working here, A'^^ = 1, A'c = 8, the renormalization group 
equation given by Arc = AGutexp[-16ir'^ /2boglut] = 9 x W^^GeV with Agut = W^^GeV, g^^t = 
47r/25.7 the unification scale and coupling |]l9|, respectively, and bo = SNc — Nf = 23 the one- 
loop beta function coefficient. It is clear that Arc 7^ Ac and so the model cannot be unified 
with the SM groups. If we insist in gauge coupling unification we would need, on top of the 
original 1 chiral -|- 1 antichiral fields, 37 extra chiral fields. If all extra fields are chiral then 
they would not contribute to the ADS superpotential in eq.(§). Of course, we think that such 
a model is not natural but it is possible (4-D string models can have different number of chiral 
and antichiral fields). 

The acoustic scale, that sets the scale of acoustic peaks, for a flat universe is given by Q 

Ia = TT^^^ (7) 

CsTls 

where To and Tis are the conformal time today and at last scattering (r = J dta~^{t),a{t) the 
scale factor) and Cs = Tj^^ Jq^" cItCs is the average sound speed before last scattering (c^^ = 
3 -|- {9/A)wb{t)/wr{t) with Wb = Ofo/i^, Wr = ^rh^ the fraction of baryon and radiation energy 
density, respectively). The acoustic m-th peak Im is then given in terms of I a and a peak 
and model dependent phase shift ip^, Im = /a("i — (fm)- It has been observed in that the 
third peak is quite insensitive to different cosmological parameters that enter in determining 
993 ~ 0.341 and so the position of the third peak is a good quantity to extract the acoustic scale 
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I A- The data from Boomerang and Maxima set the first three acoustic peaks, (the first through 
1^/2) and the acoustic scale at [|| 

/i = 213;i°, /2 = 541^20, /3 = 845ti, /3/2= 41611, ^^ = 316^^. (8) 

We have solved eqs.(Q) numerically with initial conditions J70j(Ac) = 0.3 at Ac = IGeV imposing 
ho = 0.65, Q(j,o = 0.75 and we obtain 

Ia = 307, w^o = -0.68, Weff = -0.75, (9) 

where Wejf = J da Q{a)'Wcf,{a) / J da ri(a). The energy density al last scattering (LS) is negligible 
{Q(j,{LS) = 10^^) and the average sound speed at LS is Cg = 0.52. The result is not highly 



sensitive to the initial conditions and a change in ^l^i of 50% will still be ok |18]. 

We see from eq.(|9|) that the acoustic scale I a is within the observational range given by eq.(0). 
The prediction of the first three acoustic peaks and first through {I3/2) is 

h = 223, I2 = 542, ^3 = 829, ly2 = 414 (10) 

for a baryon density wi, = 0,i,hg = 0.02 and = 1 the index of power spectrum of primordial 
density fluctuations. The peak values in eq.(|^) are consistent with the observational data in 
eq.(§) and we have used the phase shifts given in |^. 

The value of I a is sensitive to Q^f^o and even more to ho. For increasing ho (with il^o fixed) we 
find a decreasing Ia, e.g. for ho = 0.6, 0.65, 0.7, 0.8 one has Ia = 309, 307, 288, 272 respectively, 
while for an increasing energy density one has an increasing Ia, e.g. 0,^o = 0-6, 0.7, 0.75, 0.8 one 
gets Ia = 291, 300, 307, 317 with fixed ho = 0.65. Since Cs depends on Wb and an increase in 
Wh makes Cs smaller and we have therefore a slight increase in the acoustic scale, e.g. for Wb = 
.019, 0.020, 0.022, 0.026 one gets a value of U = 305, 307, 309, 314 with Q^o = 0.75, ho = 0.65 
fixed. A change in does not affect Ia but it changes the acoustic peaks through the phase 
shifts ifm slightly. 

Another relevant cosmological quantity is the density contrast on scales of Sh^^Mpc, as, which 
is constraint by the galaxies cluster abundance. For a flat universe the empirical fit of different 
authors (which converge within one a) are: Eke et al. have fig = (0.52 it 0.08)r2~°-^^+°-^^^™ 
Viana et al. report their best fit at ag = 0.56il~'''^'^(l =b 0.3) |^2|, while Steinhardt et al. [p^ ] 
have 0-8 = [(0.5 - 0.1 G) ± 0.1]n-"', @ = {n - 1) + {ho - 0.65), 7 = 0.21 - 0.22Ti; + 0.259. For 
n = l,ho = 0.65 one has as = 1.02 ± .15, 1.07 ± 0.03, 0.98 ± 0.2 for Eke, Viana and Steinhardt 
respectively. Recent analysis give slightly lower values of as = (0.46^007 )^m'''^^ for [0 (see 



also 1 25 1) and depends quite strongly on Q^n (for smaller one has a larger as). The central 
value for = 0.25 is as = 0.945 which agrees quite well with our the value obtained in our 
model as = 0.95. 

To conclude, we have shown that inverse power law potentials with n > 1 are not disfavored 
with the existing cosmological data and we have shown an explicit example with n ~ 2.57. In 
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particular, the values of w<^0 5 the acoustic scale and peaks and the density contrast ug lie within 
the observed data. The model has been derived from particle physics, using ADS superpotential, 
from a non-abelian gauge group with Nc = 8, Nj = 1. Since Nf = 1 there is only one degree 
of freedom below the condensation scale, i.e. the condensate or quintessence field (j). The 
condensation scale is at IGeV a very interesting scale which connects quintessence with the 
standard model. 

This work was supported in part by CONACYT project 32415-E and DGAPA, UNAM project 
IN-1 10200. 
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Figure 1: Variations on lead to different physical situations given by w^o and ^l^o- We have 
taken = 0.5, 0.3, .2 short-dashed, long-dashed and solid lines, respectively. The vertical line 
marks the time at Cl^o = 0.7 with /iq = 0.67. Notice that w^o increases with decreasing Cl^i. 
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